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Abstract 

The Zbb coupling determined from the Z-pole measurements at LEP/SLD shows an about 3cr 
deviation from the SM prediction, which would signal the presence of new physics in association 
with the Zbb coupling. In this work we give a comprehensive study for the full one-loop supersym- 
metric effects on the Zbb coupling in both the MSSM and the NMSSM by considering all current 
constraints which are from the precision electroweak measurements, the direct search for sparticles 
and Higgs bosons, the stability of Higgs potential, the dark matter relic density, and the muon g — 2 
measurement. We analyze the characters of each type of the corrections and search for the SUSY 
parameter regions where the corrections could be sizable. We find that the sizable corrections 
may come from the Higgs sector with light txia and large tan/3, which can reach —2% and —6% 
for pi and sin 2 6*^, respectively. However, such sizable negative corrections are just opposite to 
what needed to solve the anomaly. We also scan over the allowed parameter space and investigate 
to what extent supersymmetry can narrow the discrepancy. We find that under all current con- 
straints, the supersymmetric effects are quite restrained and cannot significantly ameliorate the 
anomaly of Zbb coupling. Compared with x 2 /dof = 9.62/2 in the SM, the MSSM and NMSSM 
can only improve it to x 2 /dof = 8.77/2 in the allowed parameter space. Therefore, if the anomaly 
of Zbb coupling is not a statistical or systematic problem, it would suggest new physics beyond the 
MSSM or NMSSM. 
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I. INTRODUCTION 



Although most of the electroweak data are consistent with the Standard Model (SM) to 
a remarkable precision, there are still some experimental results difficult to accommodate 
in the SM framework. A well known example is that the effective electroweak mixing angle 
sin 2 9 e f f determined from the leptonic asymmetry measurements is much lower than the value 
determined from the hadronic asymmetry measurements jl, 2], and the averaged value over 
all these asymmetries has a x^/dof of 11.8/5, corresponding to a probability of only 3.7% 
for the asymmetry data to be consistent with the SM hypothesis. Such a large discrepancy 
mainly stems from the two most precise determinations of sin 2 9 e f f , namely the measurement 
of Alh by SLD and the measurement of the bottom forward-backward asymmetry A h FB at 
LEP, which give values on opposite sides of the average and differ by 3.2 standard deviation. 
It is interesting to note that if such a discrepancy is attributed to experimental origin and 
thus the hadronic asymmetry measurements are not included in the global fit, then a rather 
light Higgs boson around 50 GeV is indicated from the fit {3], 4], which is in sharp contrast 
with the LEP II direct search limit of 114 GeV jlj] and results in a compatible probability 
as low as 3%. If we resort to new physics to solve this discrepancy, the new physics effects 
must significantly modify the Zbb coupling while maintain the Z-boson couplings to other 
fermions basically unchanged. In this work we focus on the Zbb coupling and scrutinize the 
supersymmetric effects. 

In our analysis we choose to parameterize the Zff interaction at Z-pole in term of the 



parameter pf and effective electroweak mixing angle sin 2 9 e ^ 



3: 



rJ /7 = (V2G,p f )^m z Y [-2Q, sin 2 9 f eff + - 75)] (1) 

This parametrization is preferred from the experimental point of view because all the mea- 
sured asymmetries are only dependent on sin 2 9 e ^ and their precise measurements can di- 
rectly determine the value of sm 1 9l^. From the combined LEP and SLD data analysis, 
the fitted values of pf and sin 2 9^^ agree well with their SM predictions for leptons and 
light quarks, but for the bottom quark their fitted values are respectively 1.059 ± 0.021 
and 0.281 ± 0.016 (with correlation coefficient 0.99), which significantly deviate from their 
SM predictions of 0.994 and 0.233 (for m t = 174 GeV and = 115 GeV) and leads to 
X 2 /dof = 9.62/2 (corresponding to a compatible probability of 0.8%). To best fit the ex- 
perimental data, pb and sin 2 9\^ should be enhanced by about 6.5% and 20%, respectively. 
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While we can envisage that the supersymmetric effects are not usually so large, we want to 
figure out to what extent supersymmetry can improve the situation. For this purpose, we 
choose two popular supersymmetric models: the minimal supersymmetric model (MSSM) 
8] and the next-to-minimal supersymmetric model (NMSSM) j^]. 

For the NMSSM effects on Zbb coupling, which have not been studied in the literature, 
we will perform the calculation to one-loop level. For the MSSM effects, which have been 
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121 ]. we will renew the study in the parametrization of p\, 



studied by many authors 
and sin 2 (the previous studies usually examined the effects on the Z-width, the ratio Rb 
and the asymmetry A b FB ). For both the MSSM and NMSSM, we will consider various current 
experimental constraints on the parameter space, which are from the precision electroweak 
measurements, the direct search for sparticles and Higgs bosons, the stability of the Higgs 
potential, the cosmic dark matter relic density, and the muon g-2 measurement. 

This paper is organized as follows. In Sec. II we introduce the general formula for the 
calculation of pf and sin 2 ^j and apply them to the MSSM and NMSSM. In Sec. Ill we 
summarize the constraints considered in this work and briefly discuss their characters. In 
Sec. IV and Sec. V we perform numerical study for the corrections to pb and sin 2 
in the MSSM and NMSSM, respectively. We will first show the characters of different 
type corrections, then we will scan the whole SUSY parameter space to investigate the 
compatibility of the supersymmetric predictions of pb and sin 2 with their experimental 
results. Finally, in Sec. VI we conclude our work with an outlook on the possibility of 
solving the Zbb anomaly. 



II. GENERAL FORMULA TO CALCULATE p f AND sin 2 9 f eff 

In the SM with the input parameters the Fermi constant Gp, the fine-structure constant 
a, Z-boson mass mz and fermion masses rrif, the electroweak mixing angle Syy = singly is 
determined at loop level by 




If I 4ira 



V2G„m 2 z 1 - Ar 



(2) 



where Ar is given by 



Ar = + ° ( 6 + l-ifvv ln(1 _ s t ) J w (f ( 3 , 



with Ti W denoting the renormalized P^-boson self-energy, 5 V and 5 b being the vertex correc- 
tion and box diagram correction to p decay p —* z/ M ez/ e , respectively. To get a more precise 
numerical result for Syy, one can iterate Eqs.(j2]) an d © a few times. 

With the sw defined above, the effective Zff coupling at Z-pole takes the following form 

BQ 

T^ Zf j = (V2G^(l - Ar)J 2 m z Y {v f - a f -f 5 + 5v f - 5a f -f 5 

-\ p z (m|) + 5Z*\ (v, - a fl5 ) - 2Q /S 2 w A«j , (4) 

where Vf = // — 2QfS 2 v and df = l( are respectively the vector and axial vector coupling 
coefficients of Zff interaction at tree level, and 5vf and 5a f are their corresponding correc- 
tions. Y! z is the derivative of the unrenormalized Z-boson self-energy T, z with respect to 
the squared momentum p 2 , and 5Z Z is the field renormalization constant of Z-boson given 
by 

5Z Z - E' (0) 2 ^ - s2]V i c w- s w f ^eS z (m|) ReZ w {m 2 w ) \ 

' 2 7 s w c w m 2 z s 2 w V m l m w / 

and Ak is given by 

A/z _ cly \ Y, z (m 2 z ) Y, w (m 2 w ) s w S 7Z (m|) + E 7 z(0) \ 
\ m| c w m| J 

In Eq.(jlj) the factor |(E' z (m|) + 5Z.f ) comes from the fact that the residue of the renor- 
malized Z propagator is different from 1, while the last term enters due to Z — 7 mixing at 
Z-pole. 

If we re-express = in Eq. (j4]) in term of p f and sin as in Eq. (JTJ , we get 

p f = l + 5p se + 5p fiV , (7) 
sin 2 9 f eff = (l + 5K se + 5nf tV )sw, (8) 

with 5n se = Ak and 

S z (0) S w (0) o s w S 7Z (0) S z (m|)-S z (0) 2 

CPse — 9 2 2 1 2 ^Z^zJ) 

5p f , v = 2^ - 2cT - 5\ 
a f 

_ a f 6v f -v f Sa f 

f ' V ~ -2Q f a f s 2 w ■ {J) 



In above equations the subscript 'se' means the contribution from the gauge boson self-energy 
which is flavor independent, and '/, v' denotes the contribution from the vertex correction 
to Zff interaction. In practice, it is convenient to express 5pf >v and 5Kf >v in term of 8g^ 
and 5g R respectively 

5p ftV = 6gi ~ ^ - 25* - 5\ 5 Kf , v = ( a f- v f) 6 9l + (af 2 + VfW R (1Q) 
a f -4Q f a f s w 

where Sg{ jR = Svf ± 5aj are the corrections to Zf L f L and Zf R f R interactions, respectively. 
From above equations one can learn that the correction to Spf tV is decided by the competition 
of 8g L and Sg R , while is mainly determined by 5g R due to (a/ + — vf) ~ 5.4. 



ings in SUSY models usually differ from 
, £ 7 z and SKf jV in the above formula 



Noting that the Feynman rules for Z-boson coup 
their corresponding rules in the SM by a minus sign 
should change sign if one uses the Feynman rules in SUSY models. The self-energies and 
the vertex corrections in SUSY models then include both the SM-particle loop contributions 
and SUSY-particle loop contributions. Since the SM-particle contributions are well known, 
in Appendix A and B we only list the one-loop expressions for the SUSY contributions. The 
only subtlety one should note is to avoid the double-counting of the Higgs contributions. 
This problem arises due to the following reason. On the one hand, the SM values of p^ 
and sin 2 are known to higher orders, and one usually incorporates such high-order 
SM effects when performing numerical calculations in SUSY models. On the other hand, 
because the SUSY Higgs sector is quite different from the SM, one cannot get the SUSY 
Higgs contributions simply by adding some additional terms to the SM Higgs contributions. 
In our calculation in SUSY models, to avoid the double-counting of the Higgs contributions, 
we first subtract the SM Hig gs c ontributions from their SM values (calculated by the codes 



TOPAZO 1161 and ZFITTER 



171]). and then we add the full one-loop contributions from the 



SUSY Higgs bosons and sparticles. 



III. CONSTRAINTS ON SUSY PARAMETERS 

Before we proceed to discuss the SUSY corrections to Zbb coupling in the MSSM and 
NMSSM, we take a look at the SUSY parameters involved in our calculations. From the 
expressions of Zff vertex correction listed in Appendix B, one can learn that the SUSY- EW 
correction depends on the masses and the mixings of top squarks, bottom squarks, charginos 
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and neutralinos, the SUSY-QCD vertex correction depends on gluino mass and the masses 
and the chiral mixing of bottom squarks, and the Higgs-mediated vertex correction depends 
on the masses and the mixings of Higgs bosons. The expressions of the gauge boson self- 
energies listed in Appendix A indicate that the SUSY correction also depends on the masses 
of sleptons and the first-two generation squarks. About these SUSY parameters, we consider 
the following constraints 

(1) Constraints from the direct search for the sparticles at LEP and Tevatron 18] 

rn^o > 41 GeV, rn^o > 62.4 GeV, m^o > 99.9 GeV, m^± > 94 GeV, 
m g > 73 GeV, m A > 94 GeV, m f > 81.9GeV, m~ q > 250 GeV, 
m~ t > 89 GeV, m~ b > 95.7 GeV, m- g > 195 GeV, 

where m^o denote the masses of the neutralinos and rrig denotes the masses for the 
first two generation squarks. 



(2) Constraint from the direct search for Higgs boson at LEP 19] . This constraint can 
limit the values of tjia, tan /3 and the masses and the chiral mixing of top squarks. In 
case of large tan (3, it can also put constraints on the masses and the mixing of bottom 
squarks. Generally speaking, this constraint re quir es the product of two top squark 
masses, m^m^, should be much larger than m 2 [20]. 

(3) Constraint from the theoretical requirements that there is no Landau pole for the 
running Yukawa couplings Y b and Y t below the GUT scale, and that the physical 
minimum of the Higgs potential with non-vanishing (H u ) and (Hd) is lower than the 
local minima with vanishing (H u ) and (Ha). 

(4) Constraints from precision electroweak observalbes such as pi ep t, sin 2 O^f}, Pc, sin 2 Q c e ^ 
and M\y- These constraints are equivalent to those from the well known e<(i = 1, 2, 3) 
parameters 23| or S, T and U parameters 24]. The measured values of these observ- 
ables are [l| 

pi ept = 1.0050 ± 0.0010, sin 2 6 l $ = 0.23153 ± 0.00016, 

p c = 1.013 ± 0.021, sin 2 6 c eff = 0.2355 ± 0.0059, M w = 80.403 ± 0.029 GeV, 
and their SM fitted values are pgf = 1.0051, sin 2 9 l e J}' SM = 0.23149, p s c M = 1.0058, 

1 eff 



sm 2 6 c eff = 0.2314 and M w = 80.36 GeV for m t = 173 GeV and m h = 111 GeV. In 
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our calculations we require the theoretical predictions to agree with the experimental 
values at 2a level. 



(5) Constraint from Rb = YiZ — > bb)/T(Z — > hadrons). The measured value of Rb is 
tC v = 0.21629 ± 0.00066 and its SM prediction is R§ M = 0.21578 for m t = 173 GeV 



ySUSY 



18|. In our analysis, we require Rf is within the 2a range of its experimental 



value. 

(6) Constraint from the relic density of cosmic dark matter, i.e. 0.0945 < VLh 2 < 0.1287 
211 ] . This constraint can rule out a broad parameter region for guagino masses M 12 , 



p parameter, tua and tan/? [22 ] . 



(7) Constraint from the muon anomalous magnetic momentum, a M . Now both the theo- 
retical prediction and the experimental measurement of a M have reached a remarkab 



precision, which show a significant deviation a e * p — a SM = (29.5 ± 8.8) x 10 10 25] 



In our analysis we require the SUSY effects to account for such difference at 2a level. 

Note that in our analysis we do not include the constraints from B physics, like b — > s'j 
281 ] and B s — B s mixing [29j, because these constraints are sensitive to squark flavor mixings 
which are irrelevant to our discussion. 

Among the constraints listed above, the constraints (4) and (5), especially the observables 
M w , pi ept , sin 2 Ogjt and Rb, are most relevant to our study of pb and sin 2 9 h e ff- Let us look 
at these constraints in more details. 

First, the precise measurements of M w , pi ept and sin 2 O^ft stringently constrain Sp se , 5n se 
and the gaugino loop contributions to Spb tV and 5nb,v The approximate forms of the SUSY 



corrections to M\y, 8p se and SK se [26( in case of heavy sparticles are given by 



5M W 


b W 


5(Ar) 


cw Ap 


M w 




[1 - Ar) 


r 2 - s 2 2 


5p se 


- Ap, 






OK se 










b w 







(11) 



where 



Ap = - _ 2 Sln9w S ^ (0) (12) 

m z m w cos^ m| 
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is the correction to the classical p parameter 6| and is only sensitive to the mass spectrum of 
the third generation squarks. Through the above relations the precisely measured M\y then 
stringently restricts Ap (of order O(10~ 4 )) and subsequently restricts Sp se and 5tz se . This 
restriction together with the precisely determined pi ept and sin 2 9 l ^Jl stringently constrains 
the magnitude of 5pi )V and 5ki >v defined in Eq.flU]) to be below O(10 -4 ). Since the gaugino 
loop effects in 5pb, v and 6nb,v are strongly correlated with 5pi >v and Ski jV (the main difference is 
caused by the mass difference between sleptons and squarks), the gaugino loop contributions 
to Spb tV and Sk^v are also suppressed, which are found to be below 5x 10~ 4 from our numerical 
calculations. 

For the constraint from the precision observable Rb, an interesting character is that it 
does not stringently constrain the magnitude of 5vb and 5cib, but it favors the relation 
Svb ~ — 1.445a&, which can be seen from the expression of the radiative correction to Rb 

lialulbi 

6Rb ~ 2 ( o 7^4 h( 3 - P 2 ) 5v t> + 2a b P 2 5a b ] oc (5v b + lM6a h ) (13) 



with p = a/ 1 — ml I m| being the velocity of bottom quark in Z decay. 

Now we turn to the constraint from the muon anomalous magnetic momentum. To get 
an intuitive understanding of this constraint, we look at a simple case of the MSSM that all 
the gaugino masses and soft-breakin g m asses in smuon sector have a common scale M. In 



this case, a^ USY is approximated by [27] 



100 GeV\ 2 



a bUbY _ 13 x 1Q -io / _ j tan/3 signfay (14) 

The gap between a^ M and a e * p then prefers a positive p, and constrains the product 
( 10 °M eV ) 2 t an P m the range [1.0,3.6] at 2cr level. So the SUSY scale can be higher for 
larger tan (3. 

In our calculations we use the code NMSSMTools [30I to generate the masses and the 
mixings for all sparticles and Higgs bosons in the framework of the NMSSM with all known 
radiative corrections included. There are two advantages in using this code. One is that all 
the masses and the mixings in the MSSM can be easily recovered if we set the parameters 
A = k ~ and A R to be negatively small. The other is that it incorporates the code 
MicrOMEGAs [3l| which calculates the relic density of cosmic dark matter. It should be 



noted that the current version of NMSSMTools only includes the constraints (1), (2), (3) 



and (6), and we extend it by including the constraints (4), (5) and (7). We note that the 
muon anomalous magnetic momentum was recently calculated in the NMSSM 32] and our 
calculations agree with theirs. 



IV. ONE-LOOP CORRECTIONS TO p b AND sin 2 9 b eff IN MSSM 

In this section we investigate pb and sin 2 to one-loop level in the MSSM. As discussed 
above, the self-energy corrections to these two observables are generally small and thus we 
mainly scrutinize the vertex corrections which include the SUSY-EW corrections, the SUSY- 
QCD corrections and the Higgs-mediated vertex corrections. We pay special attention to 
the cases where the magnitudes of the corrections are large, and show that tan (3 is crucial 
in enhancing the vertex corrections. Our analysis is organized as follows: we first investigate 
the characters of the vertex corrections to get an intuitive understanding of them, then 
by scanning over the MSSM parameter space, we study the compatibility of the MSSM 
predictions for pb and sin 2 #grr with their experimental results. 



The SM input parameters involved in our calculations are taken from [18|, which are 
a = 1./128.93, G F = 1.16637 x 1(T 5 , a s (m z ) = 0.1172, m z = 91.1876 GeV, m b (m b ) = 4.2 
GeV and m t = 172.5 GeV. 

A. Characters of vertex corrections in MSSM 

As for the SUSY-EW contribution to 8pb, v and Sk^v, the parameters involved are guagino 
masses M 1)2 , Higgsino mass /i, tan/3 = v 2 /vi with t> 12 being the vacuum expectation values 
of the Higgs fields, the soft-breaking masses Mq 3 , M Ua , M Da , and the coefficients of the 
trilinear terms A t and A b . The first four parameters enter the mass matrices of neutralinos 
and charginos, and the last seven parameters affect the masses and the chiral mixings of the 
third generation squarks [8j. 



As discussed in the preceding section, the gaugino loop contribution is small, and hence 
we only discuss the Higgsino loop contribution. The magnitude of such Higgsino loop con- 
tribution is sensitive to tan /3, the Higgsino mass p., and the masses and the chiral mixings 
of the third generation squarks. There are two characters for this contribution. One is that, 
due to the fact that the bottom Yukawa coupling Yb is proportional to 1/ cos j3, the contri- 
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bution can be potentially large in case of large tan (3 and small p. The other is that the 
contribution is moderately sensitive to the chiral mixings of the third generation squarks, 
and potentially large contribution comes from the case where the mixin g is small and the 



component of the lighter squark is dominated by the left-handed squark ll|. To illustrate 
these characters we consider three cases in the squark sector: 

(I) M s = M Qs = M U3 = M D3 = 400 GeV, A t = A b = 800 GeV; 
(II) M Qs = 200 GeV, M U3 = M D3 = 600 GeV, A t = A b = 800 GeV; 
(III) M Q3 = 600 GeV, M Uz = M Ds = 200 GeV, A t = A b = 800 GeV, 
and fix other SUSY parameters as 

M 1 = 75 GeV, M 2 = 150 GeV, m A = 500 GeV, M SU sy = 1 TeV, (15) 

where Msusy denotes the soft-breaking parameters for sleptons and the first-two generation 
squarks. Case-I corresponds to maximal chiral mixing case, Case-II is the small mixing case 
with the component of the lighter squark dominated by the left-handed squark and Case-Ill 
is also the small mixing case but with the component of the lighter squark dominated by 
the right-handed squark. 

In FigH] we show the dependence of the SUSY-EW contribution to 5p b , v and 8n b ,v on 
tan/3 in the three cases. One can see that both 5p b , v and 5K b ,v are sensitive to tan/3. As 
tan (3 increases, 5p b)V and 5n b >v get more negative contributions and, for small p, they become 
negative with sizable magnitudes. This behavior can be understood as following. As tan/3 
gets large, the bottom Yukawa coupling increases and the correction to the right-handed Zbb 
coupling 5g b R increases positively, and then Sp b)V and SK b;V get more negative contribution 
from the increasing 5g R (see Eq. tflOl) and also 5g b R in Appendix B). One also see from these 
figures that the magnitude of 5n btV is usually larger than 5p btV . The factor sin 2 9w in the 
denominator of 5K b>v (see Eq.((9]) ) can to a large extent account for this. 

Note that in these figures we only plot our results within the range of tan (3 that survives 
the constraints (1-5). The constraint (7), i.e. the muon anomalous magnetic moment, can 
in principle also limit tan (3. But this constraint relies on the mass scale of smuon, Msusy 
in Eq.l fTBl) . which p b and sin 2 are not sensitive to, so we do not apply it in plotting these 
figures. Our numerical results indicate that the muon anomalous magnetic moment allows 
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FIG. 1: SUSY-EW contributions to 5pb, v and $Kb,v with constraints (1-5). 



for a vast region of Msusy and /x where tan/3 can be as large as 60, and hence the sizable 
SUSY-EW corrections to p& and sin 2 are possible. For example, with the parameters 
in Eq.l fTB"]) . the range of tan/3 allowed by the muon (7 — 2 is tan/5 > 25 for // = 200 GeV, 
tan/5 > 33 for \i = 500 GeV, and tan/5 > 44 for p = 800 GeV. If we choose M SU sy = 0.5 
TeV, these allowed ranges are correspondingly given by 7 < tan/5 < 57, 12 < tan/5 < 71 
and tan /5 > 14. 

Next we discuss the SUSY-QCD corrections. The relevant parameters are gluino mass 
and Mq 3 , Mo 3 and Xf, = (Af, — /itan/5) which enter the mass matrix of the bottom squarks. 
From the large strength of the strong coupling, g s (mz) — 1.2 ~ 50 x Y b SM , one may naively 
postulates that the SUSY-QCD contributions to 5pb, v and Sk^^v should be much larger than 
the Higgsino loop contributions in case of rrig ~ /i and tan f3 <C 50. However, our numerical 
results show that in case of small sbottom chiral mixing the SUSY-QCD contributions to Spb iV 
and 5nb, v are negligibly small. The underlying reason is that for the SUSY-QCD corrections 

11 



g- -20 SUSY-QCD 

to : 

■25 '- m v = 200GeV 

-30 r Case-I 

-35 E-l 


5 pr 

r 
-5 - 



Case-ll 



^ - 



10 

X -15 
> 
■ef 

y -20 

to 

-25 
-30 
-35 



Case-I 



Case-ll 



Case-Ill 



20 40 60 20 40 60 20 40 60 



Case-Ill 



20 40 60 20 40 60 20 40 60 

fanp 

FIG. 2: SUSY-QCD contributions to 5pb, v and 5k\, jV with constraints (1-5). 



there is a strong cancellation between different diagrams in case of small sbottom chiral 
mixing, which can be seen from the expressions of 8g L R listed in Appendix B. It should be 
noted that such a cancellation can be alleviated for a large sbottom mixing, or equivalently, 
a large term /xtan/3 appeared in the non-diagonal elements of sbottom mass matrix (we 
checked this from numerical calculations). So the contribution may be sizable in case of 
large /itan/3, as shown in FigfJl 

Compared with the Higgsino loop corrections, the SUSY-QCD contributions in Figj2] 
exhibit a similar behavior with respect to tan/5. The difference is that the most sizable 
effects come from Case-I (maximal sbottom mixing case) with large /i, instead of Case-II 
with small u for the Hieesino loop corrections. 

n 

Finally, we consider the Higgs loop contributions to 5pb tV and 5K btV [33|. To calculate 
this part of contribution, we need to know the masses and the mixing of the Higgs bosons, 
which are determined by rriA and tan/3 at tree-level, and also by the soft-breaking masses 
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FIG. 3: Higgs loop contributions to 8p\> v and 5k}, v in Case-I with constraints (1-5). 



for the third generation squarks if the important loop correction to the Higgs boson masses 
is taken into account. As shown in FigBl the contributions exhibit a similar dependence on 
tan f3, and the significant contribution comes from the case of small tua and large tan/3. We 
checked that the results in Figj3]are not sensitive to p or Ms, and also not sensitive to the 
choice of different case ( Case-I, Case-II or Case-Ill). 

/,From the above figures one can infer that among the three types of corrections, the 
potentially largest correction comes from the Higgs loops, which can reach 2% for and 
6% for sin 2 ^//- Such large corrections reach the current experimental sensitivity since the 
current experimental measurements are p e ^ v = 1.059 ± 0.021 and sin 2 0°>™p = 0.281 ±0.016. 

Before we end this section, we would like to point out that in the large tan/3 limit the 
relic density of cosmic dark matter allows the possibility of small p or small rriA (but not 
both small). This can be seen from FigJU where we show the allowed regions in the plane of 
tan p versus p for different tua- In plotting this figure, we choose Case-I and fix other related 
parameters in Eq. (115ft . Fig H] implies that the SUSY-EW contribution and the Higgs- loop 
contribution to 5pb, v and Sk},^ cannot simultaneously reach their maximal values. 
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FIG. 4: The shaded regions are allowed by the cosmic dark matter relic density at 2a level. Other 
relevant SUSY parameters are fixed as in Case-I and in Eq. (|15|) . 

B. MSSM predictions for p b and sin 2 6 b eff 

As mentioned above, the extracted values of p^ and sin 2 9K. from combined LEP and 
SLD data analysis are respectively 1.059 ± 0.021 and 0.281 ± 0.016 with correlation coeffi- 
cient 0.99 [l|. This result is shown in FigfS] with the three ellipses corresponding to 68%, 
95.5% and 99.5% confidence level (CL), respectively. Noting that the SM predictions are 
pf M = 0.994 and sin 2 O^f = 0.233, one may infer that large positive corrections to pb and 
sin 2 6 b e ff are needed to narrow the gap between the experimental data and the SM predic- 
tion. As discussed in the preceding section, the MSSM corrections can be sizable for large 
tan /3, which, however, are negative and thus cannot narrow the gap. To figure out to what 
extent the MSSM predictions can agree with the experiment, we consider all the constraints 
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FIG. 5: The MSSM and SM predictions for pb and sin 2 6 b e fp compared with the LEP/SLD data at 
68%, 95.5% and 99.5% confidence level. The SM prediction pf M = 0.994 and sin 2 6 b £ f f M = 0.233 
is obtained with mt = 174 GeV and = 115 GeV. The MSSM predictions are from a scan (a 
sample of one million) over the parameter space. 

discussed in Sec. Ill and scan over the SUSY parameter space: 

< Mi, M 2 , M a , fx, Mq 3 , M Ua , M Ds , M A , M SUSY < 1 TeV, 

-3 TeV < A t , A b <3 TeV, 1 < tan/3 < 60, (16) 

Based on a twenty billion sample, we find the best MSSM predictions are p b = 0.9960 and 
sin 2 #gjj- = 0.2328, which give a x 2 /dof = 9.07/2 when compared with the experiment 
data. If we do not consider the dark matter constraint, the best MSSM predictions are 
p b = 0.99737 and sin 2 6 b eff = 0.2336, which give a x 2 /dof = 8.77/2. Moreover, we find 
that such a best case happens when p,,mA,vn- g ~ 1 TeV so that the three types of vertex 
corrections are suppressed. 
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V. ONE-LOOP PREDICTIONS FOR p b AND sin 2 9 b eff IN NMSSM 



A. Introduction to the NMSSM 

As a popular extension of the MSSM, the NMSSM provides an elegant solution to the 
/^-problem via introducing a singlet Higgs superfield S, which naturally develops a vacuum 
expectation value of the order of the SUSY breaking scale and gives rise to the required \x 
term. Another virtue of the NMSSM is that it can alleviate the little hierarchy problem since 
the theoretical upper bound on the SM-like Higgs boson mass is pushed up and the LEP II 
lower bound on the Higgs boson mass is relaxed due to the suppressed ZZh coupling or the 



suppressed decay h — > bb 36]. Since the NMSSM is so well motivated, its phenomenology has 



been intensively studied in recent years, such as its effects in Higgs physics [37]] , neutralino 



physics 



381 ] . B-physics 



39] as well as squark physics |40j. In the following we recapitulate 
the basics of the NMSSM with emphasis on its difference from the MSSM. 
The superpotential of the NMSSM takes the form 0, [joj 

W = teijHiH 3 d S + ^kS 3 + Key&UHi - h d e l3 Q l DE d - h e e l3 VEE 3 d (17) 

where S is the singlet Higgs superfield, and £12 = —£21 = 1- For the soft SUSY breaking 
terms, we take 

^soft = ^M 2 A a A a + ^MiA'A' + m 2 d \H d \ 2 + m 2 u \H u \ 2 + m 2 s \S\ 2 
+m 2 Q \Q\ 2 + ml\U\ 2 + m 2 D \D\ 2 + m 2 L \L\ 2 + m 2 E \E\ 2 
+(XA x e ij HiH 3 d S + h.c.) + (^kA k S s + h.c.) 

+{h u A u e lJ Q l UE 3 u - hdAoEi^bHi - h e A E e l] UEH 3 d + h.c.) (18) 

With the above configuration of the model, the fi parameter is given by fi = \(S) with 
(S) being the vacuum expectation value of S field, and the itla parameter in the MSSM 
corresponds to the combination m 2 A = -^^{Ax + ^) (see Eq. (l20l) ). So compared with 
the MSSM, the NMSSM has three additional input parameters A, k and A K . These three 
parameters should be subject to the constraints listed in Sec. Ill, and the argument that 
the NMSSM should keep perturbative up to the Planck scale requires A and k to be smaller 
than 0.7. 

The differences of the NMSSM and MSSM come from the Higgs sector and the neutralino 
sector. In the Higgs sector, now we have three CP-even and two CP-odd Higgs bosons. In the 
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basis [Re(H®), Re(Hj), Re(S)], the mass-squared matrix entries for CP-even Higgs bosons 
are 0, 0] 

•Msn = m A cos2 + m l sin 2 (3, 
M 2 S2 2 = ^4 sin 2 /5 + m| cos 2 /?, 



Ml l2 



sin 2 2/5 - sin 2(3 + ^{XA K + 4«/i), 
4/i z 2 A z 



(2AV 



ITLr 



m A ) sin (3 cos /5, 



At- 



•^5,13 = 2A/ii> sin /5 — — — m A sin 2/5 cos (3 — k/j,v cos (3, 
A^523 = 2A/iw cos /5 — sin /? sin 2/5 — /t/if sin /?, 



(19) 



and for the CP-odd Higgs bosons, their mass-squared matrix entries in the basis [A, Im(S)] 
with A = cos (3 Im(H®) + sin (3 Im(H®) are 



M 2 



.4 1 



sm 2(3 A 
3 A 2 t> 2 

-Xkv 2 sin 2/3 H ;r m i sin 2 2(3 — 3—fiA K , 

2 4/r A 



2^ 



m 2 A sin 2/5 — 3tzp,v. 



(20) 



Eqs. (ll9p and (|20|) indicate that the parameters A and Kfi affect the mixings of the doublet 
fields with the singlet field, A K only affects the squared-mass of the singlet field, and in the 
limit A, k — > 0, the NMSSM can recover the MSSM. One can also learn that in case of small 
A and k so that the mixings are small, the physical state with the singlet being the dominant 
component should couple weakly to bottom quarks and thus its loop contribution to p b and 
sin 2 6 b e ff should be small. 

The NMSSM predicts five neutralinos, and in the basis (— — i\ 2 , ip% i( s ) their mass 
matrix is given by {9], [30} 

^ Mi m z sin 6 W sin (3 —m z sin 6 W cos (3 \ 

M 2 — mz cos 9w sin (3 mz cos 9w cos (3 

-/i -\v cos (3 ■ (21) 

— Awsin/5 

V 2 > / 

This mass matrix is independent of A K , and the role of A is to introduce the mixings of i[) s 
with ip® and ip^, and k/i is to affect the mass of ip s . Quite similar to the discussion about 
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FIG. 6: The NMSSM electroweak contributions to 5pb iV and <5k& )?; , compared with the corresponding 
MSSM contributions ( thin horizontal lines), under constraints (1-5). 



the Higgs bosons, in case of small A, the correction to p^ and sin 2 should be insensitive 
to the value of k\x. 



B. NMSSM correction to p b and sin 2 9 b eff 

We first look at the SUSY-EW corrections in the NMSSM. Compared with the corre- 
sponding MSSM corrections, the NMSSM effects involve two additional parameters A and 
k. As discussed below Eq. (1211) . in case of small A, the corrections are insensitive to k (our 
numerical results verified this conclusion), and thus here we mainly study the dependence 
on A. We choose a value for k so that the allowed range of A is wide. 

In FigJH] we show the SUSY-EW contributions to 5p^ v and 5k^ v as a function of A, in 
which tan/? = 40, k = 0.4, A K = —100 GeV and other parameters are same as in FigJU 
One character of this figure is that both dpb^ and become more negative with the 
increase of A, which enlarges the gap between the theoretical values and the experimental 
data. Another character of this figure is that the contributions are less sensitive to A when 
p becomes large. This can be explained from Eq. fliZT]) which shows that the mixings between 
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FIG. 7: Same as Fig. [6J but for the Higgs loop contributions. 

ip s and the doublets (-0°, become negligiblly small for sufficiently large \i and thus reduce 
the sensitivity of the contributions to A. 

We now turn to the Higgs loop contributions to 8pb, v and 8nb, v in the NMSSM. For these 
contributions, besides and tan f3, the parameters A, k and A K are also involved. Noting 
that these contributions are more sensitive to A and k than to A K , we only study their 
dependence on A and k. 

In Figj7]we show the contributions versus A, where tan/5 = 40, k = 0.4, A K = —100 GeV 
and other parameters are same as in Figj3l This figure shows the same behavior as in FigJfH 
and the dependence on A becomes rather weak in case of large tua- 

In FigJHl we show the dependence of the contributions on k, as shown. This figure 
exhibits the similar behavior to FigJTl Compared with FigJT] and FigJHl one can learn that 
the contributions have a stronger dependence on A than on k. 

Like in FigJHl we a l so investigate the extent to which the NMSSM predictions can agree 
with the experiment by scanning over the SUSY parameter space in the region of Eq.( JT6l) 
and 

A, k < 0.7, -1 TeV < A K < 1 TeV. (22) 
Our result is shown in Figj9j Compared with FigJHl °ne can learn that the NMSSM cannot 
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improve the agreement and instead may exacerbate the agreement in a large part of the 
allowed parameter space. 

If we define a quantity F(X, k) —F(0, 0) with F denoting either <5p fe „ or 8k^ v with F(X, n) 
being the value of F in the NMSSM with arbitrary values of A and k, and F(0, 0) being 
the value of F in the MSSM limit, then by studying various cases we find this quantity is 
generally smaller than 5 x 1CT 3 , which means that in the allowed region for A and k, NMSSM 
only slightly modifies the MSSM predictions of p& and sin 2 Q\t*- 



VI. CONCLUSIONS 

The Zbb coupling determined from the Z-pole measurements at LEP/SLD deviate sig- 
nificantly from the SM prediction. In terms of p& and sin 2 O^ff, the SM prediction is about 
3cr below the experimental data. If this anomaly is not a statistical or systematic effect, 
it would signal the presence of new physics in association with the Zbb coupling. In this 
work we scrutinized the full one-loop supersymmetric effects on Zbb coupling in both the 
MSSM and the NMSSM, considering all current constraints which are from the precision 
electroweak measurements, the direct search for sparticles and Higgs bosons, the stability of 
Higgs potential, the dark matter relic density, and the muon g-2 measurement. We analyzed 
the characters of each type of the corrections and searched for the SUSY parameter regions 
where the corrections could be sizable. We found that the potentially sizable corrections 
come from the Higgs sector with light and large tan/3, which can reach —2% and —6% 
for p b and sin 2 O^ff, respectively. However, such sizable negative corrections are just opposite 
to what needed to solve the anomaly. We also scanned over the allowed parameter space and 
investigated to what extent supersymmetry can narrow the discrepancy between theoretical 
predictions and the experimental values. We found that under all current constraints, the 
supersymmetric effects are quite restrained and cannot significantly ameliorate the anomaly 
of Zbb coupling. Compared with x 2 /dof = 9.62/2 in the SM, the MSSM and NMSSM can 
only improve it to x 2 /dof = 8.77/2 in the allowed parameter space. 

In the future the GigaZ option at the proposed International Linear Collider (ILC) with 



an integrated luminosity of 30 ftr 1 is expected to produce more than 10 9 Z-bosons 41] 
and will give a more precise measurement of Zbb coupling, which will allow for a test of 
new physics models. If the anomaly of Zbb coupling persists, it would suggest new physics 
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beyond the MSSM and NMSSM. One possible form of such new physics is the model with 
additional right-handed gauge bosons which couple predominantly to the third generation 
quarks 42j. These new gauge bosons usually mix with Z and W so that the Zh^n and 
Wb^ifi couplings in the SM may be greatly changed. A careful investigation of top quark 
processes at the LHC, such as top quark decay to the polarized W boson 43j, may test this 
model in the near future. 
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APPENDIX A: GAUGE BOSON SELF-ENERGY IN NMSSM 

In the NMSSM the contributions to vector boson self-energy come from the loops medi- 
ated by the SM fermions, gauge bosons, Higgs bosons, sfermions, charginos and neutralinos, 
respetively. In the following we list the expressions for pure new physics contributions, 
namely from the loops of Higgs bosons, sfermions, charginos and neutralinos, respectively. 



We adopt the convention of 



301 ] for the SUSY parameters. 



^1) Higgs contribution: 

The NMSSM has an extended Higgs boson sector with a pair of charged Higgs bosons 
H , two CP-odd Higgs boson a, and three CP-even Higgs boson hi. The Higgs contri- 
bution to gauge boson self-energy arises from VHH, VVHH and VVH interactions 
and because we choose 't Hooft-Feynman gauge to calculate the contribution, the 
gauge boson contribution and the Higgs contribution are in general entangled. In our 
calculation, we are actually interested in the difference between the contribution from 
the NMSSM Higgs sector and that from the SM Higgs sector (see the discussion in 



the last paragraph of Sect. II). Since the SM contribution is well known \l4. Il5l |. we 
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only list the NMSSM contribution. 



Bs(p, m H +,m H +), 

B 5 {p, m H +,m H -i 



167T 2 

1 eg cos 29w 
16vr 2 2cos# 



w 



(Al) 
(A2) 



Z T ZZ (P 2 



'WW 



(P 2 



1 



g 1 



16vr 2 4 cos 2 6 



(\Sn\ 2 + \S i2 \ 2 )A{m hi ) + \Pi\ 2 A(m ai ) + A(m z ) 



-4 1 sin / gS , i2 - cosf3S il \ 2 \P' jl \ 2 B 22 (p, m a m hi 



-4| cos/3S i2 + sin (3Sn\ 2 B 22 (p,m z ,m hi ) 



+2 cos 2 26» 



-4 (m H + )-2B 22 (p, m H + , m H + ] 



+4m|| cos{3S i2 + sin/3Sii| 2 .Bo(p, m z , m ft 

i <? 2 r r 

— -^-<^ A(m^+) + (l^il 2 + |£ i2 | 2 )A(m fe J + A{m w ) 
-4 1 sinf3S i2 - cos f3Sn\ 2 B 22 (p, m H +,m ht ) 
-4 1 cosf3S i2 + sin/3S , il | 2 5 2 2(p ) m w ,m hi 



+ 



A(m ff+ ) + |i* |M(m ai ) -4|i^ 1 | 2 5 22 (p,m ff+ ,m a 



+4m w \ cosf3S i2 + sin f3Sn\ 2 Bo{p,m w ,m hi ] 



(A3) 



(A4) 



In above equations, g is the SU(2) gauge coupling, and S and P' are the rotation 



mass matrices defined in the Appendix A of [30| to diagonalize CP-even and CP-odd 
Higgs mass matrices, respectively. A and B 22 are the standard one- and two-point 
loop functions firstly defined in [34J. B 5 is related with standard loop functions by [35] 



B 5 (p, mi, m 2 ) = A(mi) + A(m 2 ) - 45 22 (p, mi, m 2 ). 



(A5) 



(2) Sfermion contribution: 

The sfermion contributions are given by 

1 9 2 



J WW 



(p 2 



16vr 2 2 



(A6) 
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-2sl?I 3f Q f R f a lR f al B 5 (p,m f ~ a ,m f ~J + s^Q 2 f B 5 (p,m f ~ a ,m f ~J k (A7) 



£?> 2 ) = I ^ 5 -C / gJflfe(p,m /al m /a ) > (A8) 

= l^^^ 1 *' ' 1 ^ " ^j^m^m^), (A9) 
where the color factor C/ is 3 for squarks and 1 for sleptons. The electric charge Qf is 
given by 2/3, —1/3, 0, —1 for u, d, I, respectively. I 3 f denotes the third component 
of the weak isospin, which is +1/2 and —1/2 for the up- and down- type sfermions, 
respectively. R is the rotation matrix to diagonalize sfermion mass matrix. 

(3) Chargino and neutralino contribution: 

For a generic interaction between a vector boson and two fermions, it contributes to 
vector boson self-energy in the form: 

+G?J^J^ (A10) 



where g L l R 3 is the coupling strength of the vector boson with left-handed or righ- 
handed fermions. The functions B 3 and B± are related with the standard two-point 



functions by [351] 



B 3 (p,m 1 ,m 2 ) = -B l (p,m 1 ,m 2 ) - B 2 i(p,m 1 ,m 2 ), 
B 4 (p,m ll m 2 ) = -mjB 1 (p,m 2 ,m 1 ) - mlB l (p,m lj m 2 ). 

For the charginos and neutralinos, the coefficients of their interactions with vector 
bosons take following forms: 

xtxjz = _9 { _y ii y* i _ \ Vi2 v; 2 + ^ sin 2 6 W ), gf^ 1 = -eSy 

COS U\y 



1.1 



9r Xj Z = c Jg w ^~ U ^ U n ~ \ U aUj2 + $ij sin 2 6 W ), g% Xj 7 = -e<%. 
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But as for the contribution from neutralino sector, one should note that, due to the 
Majorana nature of neutralinos, an addition factor | should be multiplied when using 
above formulae to get neutralino contribution to Z-boson self-energy. 



APPENDIX B: VERTEX CORRECTIONS TO Z —> ff IN NMSSM 

In this section we present the expressions of the radiative correction to Zff vertex in the 
NMSSM, namely Svf and 8a f defined in Eq.(j4]). In our calculation we neglect terms propor- 
tional to fermion mass except for f — b (bottom quark) where we keep terms proportional 
to bottom quark Yukawa coupling, Yj, ~ -^g, since those terms may be enhanced by large 
tan/?. Throughout this section all Z-boson coupling coefficients, such as Svf and Saj, are 
defined so that the common factor e/(2 sin 8w cosOw) has been extracted. 

To neatly present Svf and 5a f, it is convenient to introduce the quantities <5g{ with 
A = L,R, which denote the vertex correction to Zfxfx interaction and are related with 8vt 
and 5af by 5vf = {5g L + 8g^)/2 and 5a$ = (<5g{ — Sg{ i )/2, respectively. 5g{ is given by 14] 



ss{ = r A (mi) - 9l } %M)) - ^°/ !° s !' yS l? ) . Pi) 

bill V'Sft/ ''^Z 

where Tf x is the unrenormalized vertex correction to Zfxfx interaction, the second term on 
the RHS denotes the counter term arising from the fermion fx self-energy, and the last term 
is the counter term from the vector boson self-energy. 

Assuming the interaction between scalars fa with Z boson takes the form = 
gfifaz^p^ _|_p^^ we can wr it e down Tif x (mj) and the vertex function Tf x (q 2 ) mediated by 
a fermion if> and a scalar in a compact generic notation as 

(47r) 2 S /A (pJ) = C g gf m (b + (p f , m^m^), (B2) 

(Mfy.fs{<r) = -CA(<,T> ) gf^lgf^rn^Co 




-q\C l2 + C 23 ) - 2C 24 + ^ J (Pf, Pf, , m 0fc , ) 
gf m J gf ' V * Z W 2A ( Pf , p f , , , m*) J . (B3) 

Here C g is 4/3 for the gluino contribution (ip = gluino) and 1 for the others. The chirality 
index —A follows the rule: —L = R,—R = L. 
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If / is a lepton, the following combination of {ip, 0} contribute to the vertex: 



• Chargino correction: 

{*l>,<f>} = {X~,v} ■ 

9 L = -aVjv 9r =0; g = -1; 

<?P~ Z = 2(f/A^i + - % sin 2 M; 

g f* Z = 2(Vj 1 V? 1 + ^ i2 y/ 2 - Sij sin 2 (B4) 

• Neutralino correction: 

{^,0} = {X°,O: 

= _|^ i(A ^ 2 + tan ^ A g i) . ^"S = -yfigF^toaOwNrt 

9?' Z = -N^N*, + N j3 N* 3 ; gf® Z = N*,N l4 - N* 3 N t3 ; 

g W = (i _ 2 sin 2 e w )R l al R% - 2 sin 2 6 w R l a2 R%- (B5) 

If / is the bottom quark, the following combination of {?/>, 0} contribute to the vertex: 

• Chargino correction: 

{<M} = {rJ}: 

= (-1 + I sin 2 e w )R i al R% + i sin 2 6 w R l a2 R%- (B6) 

Note that in order to write the couplings in a neat form, we define Y t = 
m t /\/2mw sin (3, and Yb = rrib/\/2mw cos (3. Such definitions differ from their con- 
ventional definitions by a factor g. We adopt such a convention throughout our paper. 

• Neutralino correction: 

{^,0} = {X°,6}: 

v°hh* ( \/2 1 \ 

9l 3 " =9[^yR b ai(N; 2 - -tan^JVy - nft^j; 



= (1 - | sin 2 - | sin 2 9 w R b a2 R%; (B7) 
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• Gluino correction: 



{^,0} = {^}: 

gf» = -V2g s Rl i; gf« = V2g s Rl 2 ; (B8) 



• Charged Higgs contribution: 



ftM} = {*,#"}: 

ib(H-)* gm t tb(H-y gm b 
g L — — -= cotp; g R y — — -= tanp; 

V2mvK " y2mw 

^ = -(i-lsin 2 M; gf = \^e w] 

gi^-rn-z = cos29w (B9) 

• Neutral Higgs contribution: 

0} = {&,(/>, a, G )}: 

bbhi _ g m b c bbhi _ g m b Q 

gL — ~o o^* 2 ' 9r — —7, 75 bi2 'i 

Amy/ cos p zmw cos p 

Ua t _ igm b igm b 

9l ~ ~ 2m w cos (3 Pl2 ~ W 41 ^' 

yR 2m w cos/5 2m w ll H) 
g UG = W^± . g bbG° = igmb . 

^ z = (i-^in 2 M; gi z = -l^e w] 

g a*h iZ = ^ s . 2 p. 2 _ 5-.^^) = ^ Si2 sin/3 _ 5a cos/3) i^!, 

<, G °*^ Z = i(^ 2 cos/3 + ^isin/3). (BIO) 

Note that in the above formulas we did not include the contribution to 5g\ from the 
loop of {t, G~}. Such contribution alone is UV-convergent and should be attributed 
to the SM radiative effects. This situation is quite different for the neutral Higgs 
contribution where the effects of the loops of {b, G } are UV divergence and must be 
included with other neutral Higgs contribution to get an finite result. 

If / is the charm quark, the following combination of {ip, 0} contribute to the vertex: 
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• Chargino correction: 



{^,0} = {x + ,5}: 

xjxtZ = _ 2( y* Vji + ly*y j2 _ 5 .. gin 2 Qw) . 



g f^ Z = -2(Z7 il C/* 1 + \u i2 U] 2 - Sa sin 2 6 W ); 
2 o „ , _ s 2 



g r a HZ = (1 _ _ sin 2 MKi4*i - - sin 2 9 w Rl 2 R%; (Bll) 



2 

Neutralino correction: 

{^,0} = {X°,£}: 

— — ^ — K a2 ia,nU w l\ji; 



g c* a c,z = (_i + sin 2 e w )Rl 1 R% + - sin 2 9 w Rl 2 R%- (B12) 



Gluino correction: 



{^,0} = {g,c} : 



s = -V2g s Ri i; gf l = \/2g a Bfe, (B13) 



The above expressions then suffice to calculate all the Zf a f a vertex corrections Sg^. 
Summation should be taken over all non-vanishing coupling combinations, such as over the 
indices of sfermions, charginos, neutralinos, scalar Higgs and pseudo-scalar Higgs. 



APPENDIX C: NMSSM CONTRIBUTIONS TO THE ^-DECAY 

In the NMSSM the flavor-dependent correction to the decay fi — > z^ez/ e mainly comes 
from the loops mediated by gauginos, and the corrected amplitude can be written as 26] 



M = M B (l + 25 {v) + 5 {b) ) , (CI) 

where Mb is the Born amplitude, 5^ v ' is the vertex correction for either eu e W interaction 
or JiVyW interaction ( since we assume the mass degeneracy for the first two generations of 
sleptons, the two corrections are same), and 5^ denotes box diagram correction. 
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(1) Vertex corrections 

Similar to Eq. flBlj) . the correction to fif 2 W interaction can be expressed as 

ghf*w 6 M = YhWtf) _ lg^ w !^ fl (ml) + S /2 (mJ 2 )}. (C2) 

For the ev e W interaction, we have g e ^ eW = — 

(4vr) 2 £ eL (m 2 ) = \gf e " l \ 2 (B + B^m^m^.m^o) + \g% eK \ 2 (B + B x ) (m 2 , m* e , m - - ) , 
(4^) 2 S, e (m 2 e ) = \gf U ^\ 2 (B + S x ) (m 2 e , m, e , m^o ) + |^ + ^| 2 (B + B 1 )(ml e ,m gL ,m~+), 

(4ir) 2 W- = -^Tgf^ 

x ff L ' J m^om^+C + g^ J (-2C 24 + -) HPv e ,Pe, m~+, m- eL , m^o) 

x IffL X ' W ' m ^ m x7 C ' +5 « XlW/ (-2C24 + ^)| (p Ve ,Pe,m^o,mp e ,m-j) 
+2{gf e " l Ygf u ^g^ w C 2i (p Ue , p e , m, e , m^o , m %h ) . (C3) 

In the above equations, summation over % = 1 to 5 and j = 1 to 2 (x^) is implied. The 
coupling pi takes the following forms 



xJ5$w x" iX jw xjx^w xUJW g. PW p 

S-l = ~9r \ 9r =~9l I 9 ~~ ~^=, 

and for the three-point loop functions, since we take their external momentum to be zero, 
their expressions are greatly simplified: 

c„ (mi ,™ 2 ,™ 3 ) = 1 J_(i + »)Mi + .) + (i + . + t)Mi + . + 6) 



m| [ a6 (a + 6)6 

A 1 m 2 1 f-2(l + a) 2 ln(l + a 
C24 (mi, m 2 ,m 3 ) = — - 7 ln— 



4 4 /x 2 2 [ 4a6 
-36(a + 6) + 2(1 + a + 6) 2 ln(l + a + 6) 



46(a + 6) 



with a = m " 2 and 6 

m 3 

(2) Box corrections 
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The box diagram contributions to the // — > v^ev e amplitude can be expressed as 

iT = i {M{1) + M(2) + M(3) + M(4)} u el »P L v Ve u v ^P LUtl . (C4) 

Taking into account the normalization of the tree-level amplitude, —g 2 /2My V , the box dia- 
gram contributions can be written as 

s(b) = _^ky M{l y (C5) 

with each M(i) given by 

3 % 

16tt 2 M(2) = (^^)^f^(^^^)^f^ P *D 27 (m, M ,m, e ,m,-,m,o), 
16vr 2 M(3) = im^m^f*"^ 

16tt 2 M(4) = ^m i jf E5 «^H|^)l|"Ti)„K I ,m 5e) m--, m#). 

Here all the D-functions are evaluated at the zero momentum-transfer limit. Noting the fact 
that rrijx L ~ m e - t ~ ~ m^, we may write the D functions as 

_ , , 1 f -(1 + o) ln(l + o 

D (m 1 ,m 1 ,m 2 ,m 3 ) = —r 



m 



+ 



D 27 (m 1 ,m 1 ,m 2 ,m 3 ) 



ab 2 

-b(a + 6) + ((a + 6)(1 + a + b) + b) ln(l + a + b) 

b 2 {a + b) 2 
1 f (1 + a) 2 ln(l + a) 



2m 2 I lab 2 



1 + a + b)(-b(a + b) + ((a + 6)(1 + a) + b) ln(l + a + b)) 

2b 2 {a + b) 2 
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